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We study leading order perturbative corrections to the two point correlation function of 
the scalar held describing the curvature perturbation in a slow-roll inflationary background, 
paying particular attention to the contribution of entropy mode loops. We find that the 
infrared divergences are worse than in pure de Sitter space: they are power law rather than 
logarithmic. The validity of perturbation theory and thus of the effective field theory of cos- 
mological perturbations leads to stringent constraints on the coupling constants describing 
the interactions, in our model the quartic self-interaction coupling constant of the entropy 
field. If the self coupling constant is larger than some critical value which depends in partic- 
ular on the duration of the inflationary phase, then perturbation theory breaks down. Our 
analysis may have implications for the stability of de Sitter space: the quantum effects which 
lead to an instability of de Sitter space will be larger in magnitude in the presence of entropy 
fluctuations. 
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I. INTRODUCTION 



The inflationary universe scenario provides an explanation for the observed inhomogeneities in 
the distribution of galaxies and the anisotropies of the cosmic microwave background. According 
to inflationary cosmology, the primordial fluctuations arise from initial vacuum perturbations of 
the canonical variable which describes the coupled metric and matter inhomogeneities ( see [1] 
for a overview of cosmological perturbation theory). The standard computations are done at 
tree level. However, since Einstein's gravitational theory is non-linear, it is not valid to neglect 
interactions which appear at higher order in perturbation theory. Such interactions will lead to 
non-Gaussianities in the spectrum of cosmological perturbations, a topic which has received a lot of 
recent interest (see e.g. [2] for an early reference). Interactions also effect the two point functions 
via loop corrections during inflation (see e.g. [3], and they can even lead to back-reaction effects on 
the cosmological background itself, see e.g. [4, 5] for early works computing the effects of graviton 
and scalar loops, respectively). 

In order to study quantum field theory during inflation or, more ambitiously, to under- 
stand quantum gravity, it is necessary to systematically treat perturbation theory. As a non- 
renormaliziable theory, General Relativity needs an infinite number of counter-terms to get rid of 
the ultraviolet (UV) divergences (see e.g. the textbook of [6] and also the more recent articles 
[3, 7]). Since in inflationary cosmology it is assumed that UV modes originate in their vacuum 
state, the contribution of UV physics to loop corrections of correlation functions is tiny and be- 
yond the interest for observations. On the other hand, the contributions of infrared (IR) modes to 
correlation functions can be large in inflationary cosmology due to the large phase space of such 
IR modes. It is the Hubble radius H^ 1 which separates UV from IR scales. The Hubble radius is 
approximately constant (in physical coordinates) during the phase of inflationary expansion. For a 
massless field there are both UV and IR divergences. If we employ a brute cutoff renormalization 
scheme to deal with the two kinds of divergences, then the UV cutoff scale must be taken to be 
fixed in physical coordinates, while the IR cutoff is a comoving one [8]. Hence, during a period of 
accelerated expansion of space the phase space of IR modes increases. Thus, in contrast to UV 
physics, IR physics may well be lead to observationally relevant effects (see e.g. [4, 5, 9, 10] since 
the long distance modes can back-react on the semi-classical space-time or on other gravitational 
perturbation modes. There have been a lot of studies on the effects of infrared divergences on 
inflationary observables (see [11] for a recent review and e.g. [12] for some original studies). There 
is an emerging consensus that in the presence of only adiabatic fluctuations, the effect on cosmo- 
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logical observables will be small [10, 13]. However, as stressed in [14] in the context of studies 
of the effect of back-reaction of the infrared modes on the background cosmology itself, the story 
changes completely if entropy modes are present. Since matter contains many fields, it is not 
realistic to neglect the presence of these entropy modes. 

There is a second aspect of IR problems in cosmology: As pointed out in [17], perturbation 
theory may break down because of the contribution of IR modes to loop integrals. In thermal field 
theory this is called the "Linde problem" : In the thermal ensemble the soft gluons contributing as 
IR loops to Feynman diagrams cause the complete failure of perturbation theory since in fact the 
transverse gluons obtain a mass via screening. In de Sitter space-time, a scalar field with a mass 
which is tiny compared with the Hubble scale will cause the same problem as in thermal field theory 
(see e.g. [15] and also [16]). For an exactly massless field, the IR loops will be logarithmically 
divergent 1 . 

The analysis of infrared divergences in inflationary cosmology is complicated by the fact that we 
must treat both metric and matter fluctuations since they are coupled via the Einstein equations. 
This mixing must be taken into account in the computation of observables. The observable of 
particular interest is the curvature fluctuation variable £ [1, 19]. In particular, in the case of a 
massive matter field in the case of inflationary cosmology it turns out that the tilt of the curvature 
perturbations and of the gravitational wave spectrum is similar to what is obtained in the case 
of a massive scalar field in pure de Sitter space-time - thus there is a power law IR divergence in 
inflationary space-time. 

In this paper, we analyse the IR divergences arising via loop effects in the two point correlation 
function of the curvature fluctuation variable £ in the case of an inflationary space-time. We 
find a power law divergence rather than the logarithmic divergence which appears in pure de 
Sitter space-time. The IR divergences can be regulated, but the regulator is related to the initial 
condition at the beginning of the phase of inflationary expansion. Field theory in an expanding 
space-time is different from quantum field theory in Minkowski space-time: As the inflationary 
phase begins, the vacuum keeps generating UV modes deep inside the Hubble radius, modes which 
exit the Hubble radius, freeze out and are subsequently squeezed and become the IR ensemble of 
quasi-classical fluctuations which will back-react on the background [4, 5] and impact the newly- 
generated gravitational modes. Compared with the mild logarithmic divergence of pure de Sitter 
space-time, the power law IR divergences in inflationary cosmology will put stronger constraints 

1 See also [18] for an independent study of how infrared divergences in inflationary cosmology can lead to a breakdown 
of cosmological perturbation theory. 
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on the applicability or perturbative methods in inflation. 

However, it has been shown that the leading order contributions of IR divergences in inflationary 
cosmology cancel, which can be proven by applying either the in-in formalism [11] or by semi- 
classical relations [9]. Also, in the case of pure adiabatic fluctuations, IR divergences from naive 
curvature perturbations can be locally gauged away (see e.g. [20-23]), while in the presence of 
entropy perturbations the IR effects are locally measurable and can no longer be gauged away. 
Thus, the divergences from entropy modes are physically relevant [14]. Meanwhile, the entropy 
perturbations are the source of curvature perturbations and, more specifically, lead to the growth 
of curvature perturbations on super-Hubble scales. Thus, IR divergences will effect the observed 
curvature fluctuations. Hence, the fact that the curvature fluctuations are time-dependent, the 
presence of entropy perturbations, and the stronger IR divergences from loop effects will change 
the story of IR divergences in inflationary cosmology compared to what happens in the case of a 
spectator scalar field in de Sitter space. 

In this paper, we will study the power law IR divergences and the entropy perturbations in 
inflationary cosmology. In Sec. II, we will discuss the power law divergences arising in single field 
inflation models. This effect is general, and can be applied to multiple field models as well. The 
new feature in multi-field models is the presence of entropy fluctuations which we study in the 
subsequent section, We first consider massless entropy modes, and then massive ones. We derive 
non-trivial bounds on the magnitude of the self-coupling of the entropy modes which have to be 
satisfied in order for perturbation theory to remain applicable. Finally, we study the case of scalar 
matter fields which are not coupled to the inflaton. However, there is a coupling with gravity, it 
this also results in constraints on the self coupling constant if perturbation theory is to be reliable. 

II. CASE OF CONSERVED C 

In this section, we compute the leading loop corrections to the two point function of the cur- 
vature fluctuation ( for single field inflation, in which case £ is conserved on super-Hubble scales. 
As we will show a bit later, the most convenient gauge for the loop calculations is comoving gauge. 
We take into account that the Hubble parameter is not exactly constant during inflation. We show 
that this has an essential effect on IR divergences. Our discussion is based on describing space-time 
using Einstein's gravitational action, and considering matter to be composed of a set of minimally 
coupled scalar fields. 
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ILL Single Field Inflation 

Single field inflation is obtained by using the condensation of a single scalar field 4> with a non- 
trivial but very flat potential in order to drive space to expand nearly exponentially. The resulting 
space-time background is often called "(quasi)-de Sitter space", The background metric can be 
written in the form 

ds 2 = -dt 2 + a 2 (t)dx 2 =a 2 ( V )(-d V 2 + dx 2 ) , (1) 

where t is physical time, x are the comoving spatial coordinates, 77 is the conformal time, and the 
scale factor a(t) takes the form a(t) = e Ht = —7^ in the case of pure de-Sitter space. 

When introducing gravitational fluctuations, it is important to keep in mind that gravity is 
similar to a gauge field theory in that there is a redundant local symmetry, namely space-time 
diffeomorphism invariance. Out of the ten naive metric degrees of freedom, there are only six 
physical ones - two scalar, two vector and two tensor (see [1] for an in-depth discussion). The 
two tensor modes are the two polarization states of the gravitons, in an expanding universe the 
two vector modes can be neglected, and the two scalar metric and one matter degrees of freedom 
are constrained by two of the Einstein equations, leaving one physical mode. To get rid of the 
redundancy in the scalar sector, we can choose different gauges, but physics will not depend on our 
gauge choices. Uniform spatial curvature gauge is obtained by choosing the time slices such that 
the scalar metric perturbation (the component of the metric corresponding to the scalar mode) 
vanishes. In this case the scalar matter field perturbation will be non-vanishing. The spatial part 
of the metric takes the form 

gij = a 2 (Sij + hij) , (2) 

where hy is a transverse and traceless tensor representing the gravitons. 

The other convenient gauge in the case of single field inflation is the comoving gauge where 
the matter field fluctuation 5(f) is set to zero and the perturbation appears as a fluctuation of the 
spatial curvature: 

5(f> = 

9ij = a 2 [(1 + 2()5 tj + h^} . (3) 

The remaining degrees of freedom of the perturbations in the specific time slice are curvature 
perturbation £ and graviton hij. As we will argue below, for our purposes the comoving gauge has 
advantages. 



6 



At quadratic order, scalar metric fluctuations and gravitons decouple. In the comoving gauge, 
the quadratic actions of curvature perturbations and gravitons are as follows 



= l -j dtdxh [a 3 C 2 - a(d() 2 ] (4) 



and 



S 2tt = \J dtdx 3 [a*h% - a (fyhij) 2 ] , (5) 

where the subscripts s and t stand for the scalar and tensor contributions, respectively. In the 
above, we are working in Planck units, i.e. to obtain the dimensionless metric components of (3), 
we must divide the values of ( and h appearing in the above action by the Planck mass ra p \. 

At quadratic order, each Fourier mode evolves independently. The solutions of the equations of 
motion for k modes normalized to correspond to quantum vacuum perturbations on sub-Hubble 
scales are 

( k = J 1 (l + ikn)e-^ (6) 
y/2ey2k a 

and 

h±, k = -j=(l + ikr,)e- ik *. (7) 



Here e is the first slow-roll parameter, which is defined as the rate of change of the Hubble scale 
in one Hubble time: 

e- ^ ~ ^ (H) 

In de Sitter space, the Hubble scale is constant, while during inflation, the Hubble scale is 
decreasing since the inflaton is slowly rolling down its potential. The first slow roll parameter is 
required to be very small during the inflationary phase, and moreover, must be kept small for the 
duration of inflation which corresponds to at least 60 e-foldings. Hence there are other parameters 
which characterize the rate of change of e, and they are also required to be small. The second 
slow-roll parameter is 

V " ^ l</> 2 e , Q , 

n = — ~ H k = 2e . 9 

' V H<j> 2H 2 He { J 

Sometimes we require inflation to last more than 60 e-foldings, such that e is slowly changing with 
■jft ~ C(e 2 ) which is of next order of e. In this case rj ~ 2e. 
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Although the Hubble scale is not exactly constant, £ and h are good quantities to describe 
the perturbations in inflationary cosmology since they are conserved outside the Hubble radius as 
long as there are no entropy perturbations. The reason that they are conserved is the rescaling 
symmetry of the equations, as can be seen from the action which is total derivative outside the 
horizon, such that the curvature and graviton perturbations are not dynamical at all after Hubble 
radius crossing. 

Therefore, the amplitudes of Cfc and % for modes whose wavelengths are larger than the Hubble 
radius, can be described by the Hubble scale and the slow roll parameters at the time when the 
modes cross the Hubble radius: 

Ck= -=^= (l + ikr))e-^ (10) 



h ± , k = ^(l + ikrj)e-^ , (11) 



and 



where subscripts * mean that the corresponding quantities are evaluated at the Hubble crossing 
time. Although on super-Hubble scales curvature perturbations and gravitons do not depend on 
time, they are functions of k because inflation is not exactly de Sitter. 

= J?5 (12) 
H 2 

(h k h k ) ^ -p-. (13) 

The k dependence of the perturbations is characterized by spectral indices n s and n t , both of 
which are small numbers: 

"•-^"rM^'m-M^)^"-^' (I4) 

and 

n t = fcAi g(#2) =_2e (15) 

In inflationary models, the gravitational waves have a red spectrum, which means n t < 0. A blue 
spectrum is only possible if the Null Energy Condition for matter is violated or some non-standard 
sound speed is introduced 2 . Cosmic microwave (CMB) anisotropy experiments can measure n s , 
and current data indicate that n s — 1 is slightly negative. 



2 The fact that inflation leads to a red tilt of the tensor power spectrum leads to a way to falsify inflation. There are 
other cosmological scenarios such as "string gas cosmology" [24, 25] which predicts a blue spectrum of gravitational 
waves [26] 
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II. 2. Gauges 

In this subsection, we compare different gauges and the resulting different forms of the scalar 
perturbations during inflation, and we suggest that the comoving gauge is the best gauge to com- 
pute IR loops. Observables are gauge-invariant and the physics does not depend on the gauge 
choice. However, in terms of being able to reliably compute the observables, there often exists a 
choice of gauge which is particularly convenient. 

Two other gauges besides the comoving gauge have often been used in recent loop analyses of 
inflation. The first is a gauge which arises in the effective field theory of inflation in which one 
uses a Goldstone modes of gravity perturbations (see e.g. ir [27]). However, for dealing with IR 
loops this gauge has problems. The main idea of ir effective field theory is to apply the Goldstone 
equivalence principle to the inflationary perturbation theory. At high energy, the mixing between 
Goldstone mode and metric perturbation modes is not important, and then the scalar perturbation 
can be treated as a Goldstone boson it. For single field inflation models, this effective field theory 
breaks down for modes whose energy corresponds to E < e l l 2 H, it. But a large contribution to IR 
loops comes from modes in this range. 

A second frequently used gauge is the uniform density gauge C = 0. In this case, the fluctuations 
are carried by the scalar field perturbation 5(p which is non-zero. It is well known that 5(p is not 
conserved outside horizon. To calculate the tree level two point function of curvature perturbation, 
it is convenient to make use of the fact that £ is conserved after crossing the Hubble horizon. Thus, 
it can be calculated at the Hubble crossing time. With the relationship matter fluctuations S(p and 
metric fluctuations ( at this time (valid in any gauge), we can obtain the curvature perturbations 
while working in the uniform density gauge. However, in this gauge the IR loop calculation needs 
to take into account the time evolution of the integrands on super-Hubble scales, and thus involves 
an integration over the whole time of inflation. In other words, if we use 5cj), we must take into 
account the time dependence of 5<p after Hubble crossing time. 

On the other hand, working in comoving gauge in the single field model of inflation, the curvature 
perturbation Q is conserved outside of the Hubble radius. The time dependence of correlation 
functions of the fluctuation variables in this gauge is trivial. When we perform the time integration 
in loop integrals, the IR modes of ( are time-independent and can be taken out of the time 
integral. Thus, the computations are technically easier. Note, however, that the small tilt from 
the k-dependence of C is essential for the IR loop diagram integration. In the UV, the tilt of the 
spectrum is not crucial, since deep inside the Hubble radius the WKB method can be used [28]. 
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To summarize, we conclude that the comoving gauge is the best one for computing IR loops in 
single field inflation models. In multi-field models of inflation and in the presence of non-vanishing 
entropy perturbations, ( is no longer conserved. The IR divergences of inflation becomes more 
interesting as we will discuss it in the next section. 



II. 3. Interactions 



After choosing the gauge, the next step is to gather the interaction terms for loop calculation. 
In comoving gauge, there is no interactive ghost field and thus we need not consider any ghost 
terms. In single field inflation with sound speed c s = 1, we use e as second expansion parameter for 
the gravitational action. The first expansion parameter is the amplitude of the linear fluctuation 
C, and for each order in ( we can expand the corresponding terms in the interaction Lagrangian 
in terms of e. Here we just keep the leading terms in e for every kind of interaction. In addition, 
there are terms proportional to jjr which can be dealt with by redefinition of £ or hij. These terms 
are neglected here because they have nothing to do with the leading IR divergences. 

In this section, we will list all interaction terms that are relevant to one loop corrections to 
two-point function, and we also write down the higher interaction terms since they can be used 
to estimate the order of magnitude of the two loop corrections to two-point functions, or the loop 
corrections to n-point correlation function. In the second part of the section, we will give the 
method to estimate the result of the loop corrections. In the following, we will use conformal time 
because it is convenient in order to estimate the order of the loop corrections. 

The second order action is 0(e) and takes the form [1, 29, 30] 

S 2 = j d7?dx 3 ea 2 (V 2 - (dC) 2 ) (16) 
The cubic interaction term was first calculated in [31]. It is 0(e 2 ) and is given by 

S 3 = J d V dx 3 e 2 a 2 [C 2 C - (d() 2 ( - 2('d iX d i (\ , (17) 

where 

d\ = C ■ (18) 

The interaction Hamiltonian is written as 



hP = -c 3 . 



(19) 
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The leading fourth-order term is of the same order as the cubic interaction 0(e 2 ) (it was first 
computed in [32]) 

S 4 = J dtdx 3 e 2 a 2 {-Pjd 2 ^ - ^~ 2 S [C 2 + (d() 2 } ~ 3(cr 2 S) 2 - , (20) 

where 

I 9 4 ft = grs (Q.x rs _ 5(rSs)i ) (21) 

and 

z rs = d r ('d s ( + ('d rs (. (22) 

We have also used £ = 5 rs T, rs . If keep the leading order, the interaction Hamiltonian is also the 
minus Lagrangian, 

H[ I] = -U (23) 

The fifth order and sixth order terms are more and more complicated, and we just mention the 
order in e of the leading terms in these interaction terms, and we omit partial derivatives on £ in 
the expressions (this expansion was first discussed in [33]) 

5 5 ~ J dtdx 3 ae 3 ( 5 (24) 

and 

5 6 ~ J dtdx 3 ae 3 ( 6 . (25) 

For the loop corrections, S3 is subdominant compared to S4. The one loop diagram with S3 
interactions has one more vertex than the one with S4. It is thus suppressed by e 2 , but it has one 
more propagator, which gives an enhancement of l/e. Therefore, S3 is subdominant. 

To compute loop diagrams we must identify the vertex factors and propagators. Concerning 
the vertex factors, 

e 2 

S4 vertex >->■ — — ■ (26) 
H 

and 

e 2 

S 3 vertex i-> . (27) 
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Every internal and external curvature C propagator gives 

H 2 

( propagator i-> — , (28) 

and the graviton propagator (which will be relevant since there are interactions between scalar 
metric fluctuations and gravitons at higher than quadratic order) is 

h propagator i-> H 2 , (29) 

which comes from the action for free gravitons which is 

S h2 = \j d V dx*a 2 (h^ - fyhijdthij) . (30) 

Now let us consider the individual interaction terms. The cubic interaction term between two 
scalars and a graviton is 

S^h = f dvdx^hijdiCdjC . (31) 

The cubic interaction term between one scalar and two graviton is subdominant in e since it 
has the form 

%2 = J dr&c^HcPti^dr*? . (32) 

The quartic interaction term between two scalars and two gravitons is 

S C 2h2 = J dridxha 2 (—huhijdiCaA ■ (33) 

Counting the order of e is a easy way to know whether a particular interaction is dominant or 
not. Let us consider the leading contributions to the two point function of (: 




lie \/e l/e „ x (jr divergence) (34) 




H 4 

l/e l/e ~ x (IR divergence) (35) 
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1/e 




1/e 1/e 1/e ^ #4 x ( IR divergence) (36) 



H 4 

~ x (IR divergence) (37) 



~ H 4 x (IR divergence) (38) 

3 where the wavy lines stand for gravitons and the solid lines are curvature propagators. The 
dominant interaction terms at one loop are £4, S^h and S^M- This counting method can be easily 
generalized to n— point correlation functions and to higher loop corrections. The IR divergences 
are logarithmically divergent if H is constant. 

II. 4. In-in Calculation of Loop Corrections 

In Sec. II. 3, we have obtained the interactions and discussed the method to estimate the order 
of loop diagrams. In this section, we use the "in-in" formalism to perform the actual calculation of 
the loop corrections to the correlation function via evaluating the leading order loop diagrams. As 
is well known, in cosmology there is no easily defined "out" state, so "in-in" formalism is applied 
in order to calculate the observables, e.g. the correlation functions. 

The "in-in" formalism [34, 35] is similar to the real-time formalism in thermal field theory. It 
makes use of the interaction representation of quantum field theory. In this formalism, the vacuum 
expectation value of an operator O is written as 

{n\0(t)\n} = (<d\Te i $-™ Hl{t ' )dt '0 I {t)Te- i S t -™ Hl{t ' )dt '\Q) (39) 




3 In fact, because of the derivative terms in the vertex factors, there are no IR divergences in Eq.(37). 
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where \Q) is the vacuum of the interacting theory and |0) is the free field vacuum. The right hand 
side of the equation is calculated in the interaction picture. T and T denote time-ordered and 
anti-time-ordered products. 

We are considering the two point function of the curvature perturbation (. There are two kinds 
of quantum corrections contributing to this Green's function, one of which is from curvature loops 
and the other from graviton loops. The quantum correction from eq (34) is denoted as Gp~(r)), and 
the one from eq. (35) is denoted as G^rj). Using the Feynman rules discussed in the previous 
subsection they are given by 

Gj'fa) = 26 £ d m a 2 ( m )lm [c^KfM] J ^ 2p 2 S m 2 e\h q ( m )\ 2 (40) 

and 

g "^-^iH><- (4l) 

where 6 is taken to be the angle between p and q. 

If Hk is not k dependent, the final result will be logarithmically divergent. However, in in- 
flationary cosmology we must take into account the small k dependence and study how this will 
impact on the physics. Focusing first on G^ 1 we find 



>• ' Aep 3 16vr 2 J Air q ini ""' 
I 2 

initial 



nt 



H 2 1 SI! 2 ( Am ' ; ' 



4ep 3 167r 2 n t \k. 

R2 1 4H initial 

4ep 3 IQ-k 2 e 



^initial 



(42) 



where in the final step we have identified the IR cutoff length scale with the value of the Hubble 
radius at the beginning of the inflationary phase. Taking the IR cutoff momentum to be smaller 
would make our result probe the state of modes which are super-Hubble throughout the cosmolog- 
ical evolution. We have also made use of the relation (15) between nt and e. 

With the same logic and method, the term Gp(rj) can be calculated as follows: 

Gfa) = -16e 2 Re f d m a 2 ( m ) P d m a 2 ( m ) f d 3 q p 4 sin 4 6> 

-kpWPCp^OCfe)^^)^^)^^)^^)], (43) 
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which leads to 



G 1 (n) ~ 



H 2 1 ( ^Hf niUal ^Initial 



n t n s -l 



_ H \ ( r initial , initial \ ( AA\ 

~ 4ep 3 16vr 2 \ e 3e - v J ' 1 ' 

The one loop corrections from S4 do not have any IR divergences, because the terms associated 
with loops are always have a derivative acting on them. Hence, there are no integrals such as 
f dq 3 /q 3 which result, and hence there is no IR divergence. 

Notice that Gp and G^ 1 have different signs, and thus there is a possibility that the terms might 
cancel. In fact, if we consider a theory in which e is very slowly varying on a time scale of H^ 1 
compared to its actual value, which often happens in models which admit eternal inflation or in 
which inflation simply lasts a lot longer than the minimal 60 e-foldings which it needs to last in 
order to be successful at resolving the flatness and horizon problems of Standard Cosmology, then 

77 ~ 2e (45) 

and the leading order IR divergences is cancelled. On the other hand, if is of the same order 
as e, then the analysis is more complicated, and we need to take into account the detailed time 
dependence of e. 

In models in which inflation only lasts 60 e-folding the leading IR divergences do not cancel. 
However, the magnitude of the loop corrections is small because the phase space of IR modes is 
small. 



III. ENTROPY MODES 



The restriction to single field models of inflation is not very physical. We know that matter in 
the real world is described by a large number of fields. The inflaton is one additional field which 
must be added to the theory. The inflaton must also couple to other fields in order to provide 
a successful exit from the inflationary phase. The presence of other matter fields leads to the 
presence of entropy modes. Thus, when computing quantum effects on cosmological observables it 
is essential to take into account the effects of loops involving entropy modes. 

Whereas is some circumstances adiabatic loops lead to effects which can be gauged away (as 
we have seen in the previous section), effects of entropy loops cannot (see e.g. [14]). Thus, if the 
cosmological fluctuations can be influenced by IR effects from entropy loops, these effects cannot 
be gauged away. They will contribute to the effective field theory of inflationary perturbations and 
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can hence influence cosmological observations. In this section, we will study the IR divergences 
to curvature perturbation via entropy loops. We will see that these effects can be very important 
and that these loops effects can lead to a breakdown of perturbation theory if the self-coupling 
constant of the entropy modes is larger than some critical value which is much smaller than 1. 
The magnitude of the infrared effects of entropy modes depends on the mass of the entropy modes 
and on the initial conditions. We separate this section into three parts, one of which deals with 
massless entropy modes, and the second subsection considers entropy modes with a small mass. 
These two cases will give different IR divergences. And in the final part we will deal with a more 
general case in which we take the Higgs field to be the inflaton as a example. 

Several groups have studied entropy loops [3, 11] by adding a (approximately) massless scalar 
field a to the matter sector which minimally couples to gravity. The Lagrangian for a is 

L = -\d^a. (46) 

The correlation function does not pick up important IR contributions from the kinetic term, 
since the derivatives will introduce an extra power of momentum into the loop integrals which make 
the integrals more UV rather than IR divergent. Considering the simplest interaction between the 
inflaton and the entropy field, 

C = -\d^o -g<f> 2 a 2 . (47) 

it appears that there are sizeable IR divergences. However, unless g is very small, the entropy 
field a will acquire a large mass from the expectation value of the background inflaton field. The 
effective IR cutoff for entropy loops is thus the larger of either the Hubble scale (the a mass induced 
by gravitational effects) or the induced mass of the entropy mode. Thus, the IR contribution of 
entropy modes is tiny. 

However, the conclusion that the entropy IR modes contribute little to the two point correlation 
function of £ has an important loophole, which relies on the fact that the curvature perturbation 
is no longer constant outside Hubble radius. It is well-known that entropy fluctuations source 
a growing curvature perturbation. Also, the entropy self-interaction term will contribute to the 
curvature perturbations via the gravitational coupling. In order to study the effects of entropy 
perturbations, there are two ways to proceed. The first is to solve for the time dependent curvature 
perturbations from the equation of motion for £ in the presence of entropy sources and then to 
perform the loop calculation using the in-in formalism. The second approach (the one we now 
adopt), is to solve the free equation of motion for £, thus obtaining a time-independent solution, 
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and then treat the source term as an interaction term. We then use the in-in formalism to obtain 
the corrections from IR entropy modes 4 . 

III.l. Massless Entropy Modes 

It is a general result that there is mixing between isocurvature and curvature modes. One 
simple example in which it is easy to see this mixing is to consider a model with a complex scalar 
field with a "Mexican hat" potential. The adiabatic mode corresponds to the scalar field moving 
in angular direction in the vacuum manifold, and the isocurvature direction is along the radius. 
The action describing the system is as follows [38] 

S = j dr/d 3 ^ [a 2 e' 2 {R + af - a 4 Acr 4 ) , (48) 

where R represents the radius of the circle in field space, and 9 is the angle of the scalar field 
configuration. The field 9 represents the curvature perturbation. The self interaction of the entropy 
modes can have various form; here we use a Act 4 interaction without loss of generality. If the 
interaction does not contain more than two time-like or spacial derivatives, then entropy loops lead 
to a IR divergence. 

The mixing term is treated as an interaction. As we argued in Section 2, it is convenient to 
work in the gauge given by 56 = 0. To obtain the interaction Hamiltonian in this gauge, we start 
with the perturbative expansion of (48) and then (valid to first order approximation) transform 
the interaction term to the comoving gauge by the transformation 

RS9 -> a-HH) . (49) 

The interaction Hamiltonian then becomes 5 

Hi = -2a 2 6' 59'Ra -»• 4a 3 ^(V (50) 

R 

where the second arrow indicates the result after performing the gauge transformation. If the 
entropy modes are massless, we have the exact solution of the free entropy modes which start on 
sub-Hubble scales in their vacuum state: 



(1 + ikr))(T ikri . (51) 



4 The corresponding result is equivalent to the leading order in expansion with respect to the coupling between 
entropic and adiabatic modes of the result in the first approach. This perturbative approach is employed in [36] in 
evaluating the cross-correlation between adiabatic and entropic mode and in [37] and [38] in evaluating the transfer 

of non-Gaussianities from the entropic to the adiabatic sector. 

5 In multi-field models, in general there are two types of couplings between adiabatic and entropic modes (see e.g. 

[39]): ~ CC'°" an< l ~ %£C°" with the common coefficient £ denoting the tuning rate in field space (as 9 in eq.(48)). 
In this work, in order to see the qualitative effect of entropic mode, we only take into account the contribution 
from the first type. 
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Note that these mode functions freeze out and become constant in time on super-Hubble scales. 

First, the time dependent part of tree level contribution to the two point function 6 is captured 
by the two point mixing interaction term (50) and yields 



G tree = -2Re( J Hi J tf/CC) + ( J #/CC J H T ) 



( 27 r) 3 <5 3 (k - k') (^L 

v rm 

dma 3 {v2) / d^a 3 (^)<T fc (»;i)Cfc(m)fffc(»»)Cfc(%)Cfc(»;)Cfc(^) 

-oo •J — oo 

+2 f dWfai) f d % a 3 ( % )a fc ( m )CUm)Cfc(»7)a(»7)^(%)Cr(%) 

J —oo J —oo 

= (2 7 r) 3 ^(k-k')^^x[l g(-M] 2 (52) 
The tree level correction from the entropy modes is given by the following diagram and yields 



RH 3 RH 3 



H 2 /z H 2 H 2 /e „ -LEl x [log(-fe»7)] 2 . (53) 

Here, the solid lines denote curvature correlations (((} ~ and the dotted lines are entropy 
modes with propagator {a a) ~ H 2 . The interaction vertex is proportional to In this estimate, 
the extra H 3 in the denominators of the interaction terms arise from the H-dependence of the scale 
factor a 3 = — H s 5 arising in the interactions. The entropy perturbations source the curvature 
perturbation via 

C ~ 8s , (54) 

and this explains why the tree level result above contains the two factors of log (—krj) ~ t. 

So far, we have neglected the mass and self-interactions of the entropy field. For the moment, 
we will assume that the mass of the entropy modes is small and compute the effects of the self- 
interactions. To be specific, we choose a a 4 interaction 

%j = a 4 \a 4 (55) 

and will go on to study the resulting entropy loop correction to the curvature two point function. 



6 See [36] for a related calculation of the corrections to the adiabatic power spectrum due to the coupling between 
adiabatic and entropic modes. 
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It is given by the following Feynman diagram 



' H 2 \ 
/ \ 

! I 

^ 4 ' 

S X/H / 



H 2 ~ x Alog(— krj) x (IR divergence) 

~ ^7 x [log (_A;7?)]2 x ^r^f^ x log (-M • (56) 

As in Section II. 4 (Eq. (42)), the power-law infrared divergence leads to the factor of Hf nitial /e, 
where the IR cutoff length scale is set to be the Hubble radius at the beginning of inflation. The 
physical reason for having the extra term log(— krj) is that the entropy mode also sources itself and 

H 2 

the effect increases with time. The IR divergence is expresed via the term ' n ^' a ' where the IR 
cutoff length scale is (as in Section 2) set to be the Hubble radius at the beginning of inflation. 
This divergence is similar in structure to what we found in Eq. 42. However, in this case there 
can be no cancellation in the limit of very slow-roll inflation. Note that the logarithm divergences 
of pure de Sitter space has changed to a power law divergence. Summarizing, the entropy loop 
correction to the curvature two point function is given by 

G loop = 2Im( J Hj j Hj J HjC() - 2Im( J Hj J Hj(( J Hj) 

= (2.)^ ( k _ k0 \_^lL_ x [log ( _ M] 3 . (57) 

Comparing the tree level result with the one loop contribution we see that perturbation theory 
breaks down when 

A I^2^(^f^) 2 x[ lo g(-^)] > !■ ( 58 ) 

Using dynamical renormalization group techniques [40], the temporal logarithm divergences 
[log(— kr])] can be resummed. In fact, one way to interpret this result is that interaction gen- 

/ \ 2 

erate a mass of order 4 ' 



H 

After doing this, the condition for the breakdown of perturbation theory changed slightly and 
becomes 

A i u^ 2 >L (59) 

which comes from the subleading term in the loop corrections and is still sensitive to the initial 
condition. This constraint on A is severe, in particular in models with a very long period of inflation 
when the ratio of H values is large. 
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III.2. Massive Entropy Modes 

Massive entropy modes are not constant outside the Hubble radius (even at level of free field 
theory) . The free field equation of motion is 

4'-^+(W^)a fc = 0. (60) 

After choosing the Bunch-Davis vacuum at past infinity the solutions are 

v k = -ie* v+1 M^H(-r,)lT$)(-kn) for ^<\, (61) 
where H v indicates a Hankel function and 



/9 m 2 

"-Vi-55- (62) 

If — krf <C 1 (super-Hubble scales), the values of the modes are slowly decreasing in time: 

o k ^ -e^ + 5)i^r»4- (-k7])- u+ ^ (63) 
V 71 " k2 

In the case of m <C H we have v ~ | — ^ ^ , and hence 

2 2 

(-^)- v+ § ~ (a(ri)H)-^ ~ e -5ir( fr *+ lllf ( 64 ) 



Thus the entropy modes will have decayed away on a time scale of i2"t ~ 3^-. Assuming that 
m ~ 10~ 2 i?, then the decay will have occurred after 10 4 e- foldings. Thus, the effective IR cutoff 
for the massive entropy modes is 

k IR ~ #e~ 3 ^ . (65) 

If inflation does not last as long as 3^ e-foldings, the IR cutoff is the same as for massless entropy 
modes when it is set by the initial value of the Hubble radius, as we have previously discussed. 

We now first calculate the tree level correction to the curvature correlation function. The 
calculation is similar to the analysis of Eq. (52) in the case of massless entropy. We will just write 
down the leading order terms: 



G tree = (2tt) 3 <5 3 (k - k') ^ 



R 



x 



4"" 1 (cos(2™) - cos(4™)) T{vf_ ( _ fcf/) -2„ + i 



4e 2 (v- |) 2 fc 3 

(66) 



2n 2 k 3 e 2 (2v - 1) 

The first term in the bracket is the same as Eq. (52), with the log 2 (— krj) divergence being recoved 

3 
2- 



in the limit of v — > §, and the second term is zero in this limit 



20 



The one loop result is complicated, and we here only show the order of magnitude of the result 
which can be derived from the in-in formalism similar to in Eq. (57), or can be estimated as in 
Eq. (56). The result is 

1 ( H mttial \ 2 (-kr,)<-^) 



Gloop ~ Gtree X A- ^ — J - 3 (67) 

where x is a non-negative integer of the order one. Comparing with the result for massless entropy 
modes, we see that the logarithm divergence in time becomes a power law divergence. 

Comparing the magnitudes of the tree level and one loop terms again leads to a condition on A 
in order that perturbation theory remain valid. The most conservative constraint on the coupling 
constant A, independent whether the entropy mode is massless or massive, is 

A < ee^lj (68) 

where N is the total e-folding number of the inflationary period, and N' denotes how many e- 
foldings before the end of inflation that the k mode corresponding to the external line of the 
curvature perturbation exits the Hubble radius. Choosing the fiducial values A^ ~ 60, e ~ 0.1 and 
N' ~ 0(1 -7- 10), then 

A < 10~ 7 -r- 10" 6 (69) 

However, in large field models of inflation the value of N is much larger than the minimal value 
chosen to obtain the above estimate, and then the constraint (68) will be much more stringent 
than the result (69). 



III. 3. Two Loops and Higgs Field 

In this subsection, we will discuss a more general situation which is expected to occur during 
inflation. Besides the inflaton, we expect physics to contain many other scalar fields, such as the 
Higgs field(s), moduli fields, and superpartners of fermions. There is the possibility that these 
scalars do not directly interact with the inflaton but only interact with curvature perturbations 
via gravity. If the masses of these scalar fields are light compared to the Hubble scale (which we 
expect to be true in the case of the Higgs field(s)), then the modes of these scalar field will decay 
at 3^ e-foldings after Hubble radius crossing, as discussed above. For example, the decay of the 
Higgs field with mass O(100)GeV begins at ^ ~ O(10 25 ) e-foldings after the modes leave the 
Hubble horizon. Therefore, for fields with a small mass, the IR cutoff will usually be given by the 
initial Hubble radius at the beginning of the inflationary phase. 
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The interactions terms for the scalars come from self-interactions 

Sj = J d V d 3 xNVh (-Aa 4 ) = J dr/dWe 3 ^'/^) (_Aa 4 ) 

,-/d,d3, a 4 gc 2 + ^ + _i_^) Acj4 (70 ) 

Here we neglect the interactions with the graviton, and only consider the interaction with the 
curvature perturbation £, since the two point correction function of curvature perturbations are 
considered. This kind of interaction can lead to two types of loop correction to the two point 
function of £• 

We can estimate the magnitude of the loop corrections by counting powers in H and e and 
making use of the power-law IR divergences which depend on the initial value of the Hubble scale, 
as in previous subsection. The Feynman diagram and the scaling of its contribution to the two 
point function is 

/ H 2 \ 
/ \ 
I I 
\ ' 

\ A/tf 4 / 



H 2 /e / \ H 2 /e 

I \ 

I I 

\ ' 



X H 2 / ,# 4 



A— =- x (IR divergence) 



e 2 

The magnitude of the loop correction can be calculated by in-in formalism and yields 

G 2 ioop = -2Im( J HjCQ 

= (2.)3*3 (k _ k0 x A . (72) 

We can now ask what constraint on the coupling constant results. The answer depends on 
how long inflation lasts and how large H is. If we assume that there was eternal inflation in the 
early universe, and that thus the density fluctuations initially were very large (on scales which are 
obviously completely irrelevant to current observations), i.e. 



2 H 2 

initial ^ 1 fT*^ 

pj "2e(2^M|- i ' {16) 
then we can derive the following constraint on the coupling constant A: If we require that the two 
loop correction is not larger than the tree level value, then 

A < 1(T 10 x e 2 < 1(T 12 , (74) 
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where 10 is from the thermal fluctuations of CMB. Thus, if we want to trust the effective field 
theory of eternal inflation, all scalar fields with the small mass must have a self-interaction coupling 
coonstant which is smaller than 1CT 12 . 

IV. DISCUSSION 

In this paper, we have computed one loop corrections to the two-point correlation function of 
the scalar field ( which characterizes the curvature perturbations. We emphasized that the In- 
dependence of the Hubble scale during inflation will change the IR logarithmic divergences of pure 
de Sitter space into power law divergences. The power divergence implies that the IR divergences 
are sensitive to the initial value of the Hubble parameter, the natural IR cutoff scale in the context 
of inflationary cosmology. 

We argued that the best choice of the gauge to perform the computation of IR looops in single 
field inflation models is the comoving gauge, since in it the perturbation variable (namely Q is not 
time-dependent on super-Hubble scales. 

We showed how to estimate the IR divergences in one loop corrections to the two point cor- 
relation function of (, and we also used the in-in formalism to calculate the one loop diagram, 

H 2 

both of which are consistent. The correction is of the order O ( in * tial ) . In single field models of 
inflation, then if the inflationary phase lasts long, the requirement that the slow roll condition on 
e is preserved in time will lead to some cancellations of the loop corrections, and the semi-classical 
calculation breaks down only when Hf nitial ~ 1. The cancellation of the leading IR divergences in 
single field inflation models is consistent with the semi-classical relation obtained in [9]. 

However, in multiple field inflation models, the curvature perturbation ( is not conserved on 
super-Hubble scales. The entropy modes source the curvature perturbation and also self-interact. 
Neither the kinetic term of the entropy field nor the interaction terms coupling the entropy field 
to the inflaton lead to an IR divergence, but the self-interaction terms do lead to IR divergences. 
Massless entropy modes will lead to logarithm divergences in conformal time, while the massive 
(but low mass) ones will lead to power law divergences in time. Estimating the ratio of the one 
loop to tree level results leads to constraints on the self-coupling of the entropy modes which have 
to be satisfied in order to avoid the breakdown of perturbation theory. In the two fields inflation 
models, for a Act 4 interaction term (where a is the entropy field), the constraint takes the form 
A < 10~ 7 -T- 10~ 6 if we require that the effective field theory for inflation remain valid for more 
than 60 e-folding. Because this constraint is on the self-coupling, we believe that it is rather 
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generic. If inflation lasts longer than the minimal duration it must last in order to address the 
cosmological problems of the Standard Big Bang model, then the constraints on A become much 
more restrictive. In more general cases, any scalar (for example, the Higgs field) which does not 
couple directly to the inflaton field, couples to gravity. The leading order correction which such a 
field has on the curvature perturbation is from a two-loop diagram, and the requirement that the 
two-loop correction is not larger than the tree level result leads to a constraint of the self-coupling 
of the scalar field of the form A < 10~ 12 if the inflaton field starts in the region of eternal inflation. 

Our analysis has implications for the recent debate about the stability of de Sitter space. There 
is one camp [41] arguing for the stability of de Sitter space even in the presence of fluctuations 
based on performing the analysis in Euclidean de Sitter space. A second camp [42] argues that in 
a cosmological context a de Sitter phase is unstable because of large IR back-reaction effects. The 
renormalization prescription breaks the de Sitter invariance as can be seen either by imposing an 
initial space-like Cauchy surface (e.g. the initial time of the inflationary phase), or by setting the 
IR cutoff to be fixed in comoving coordinates (see e.g. [43]) for a recent comparison of these two 
camps). We have shown that in the presence of entropy fluctuations, the effects which lead the 
second camp to argue that de Sitter space is unstable will be much larger in magnitude. Thus, our 
analysis may imply that entropy modes lead to an enhanced decay rate of de Sitter space. 
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